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Abstract. We report on extensive V, R, and i band 
photometry of the black-hole candidate Nova Sco 1994 
(GRO J1655-40) obtained during March 1996, when the 
source was close to its quiescent pre-outburst optical 
brightness (V=17.3). From our observations and data 
taken from the literature we derive a refined ephcmeris 
for inferior conjunction of the secondary star: HJD 
244 9838.4198(52) + 2.62168(14) x N. We have modeled 
the V, R, and i band light curves in terms of a Roche 
lobe filling secondary and an accretion disk around the 
compact star, the latter described as a flat cylinder with 
a radial temperature distribution. From the shape of the 
light curve we constrain the binary inclination and mass 
of the secondary star to lie in the ranges 63? 7-70? 7 and 
1.60-3.10 Mq, respectively. This limits the mass of the 
black-hole to the range 6.29-7.60 M Q . The mass range we 
obtained for the secondary star is supported by the results 
of stellar evolution calculations. 

Key words: (Stars:) binaries: close, individual: Nova Sco 
1994 (GRO J1655-40) - X-rays: stars 



1. Introduction 

Soft X-ray transients (SXTs) are low-mass X-ray binaries 
consisting of a neutron star or black-hole primary, which 
undergo unstable accretion from a late-type companion. 
Disk instabilities ( Van Paradijs 1996 ) cause brief but vio- 
lent outbursts, typically lasting weeks to months, during 
which their X-ray luminosity abruptly increases by several 
orders of magnitude, to ~ 10 37 — 5 x 10 38 erg s _1 . Dur- 
ing such an outburst, the optical flux becomes dominated 



by the contribution of the accretion disk and is enhanced 
by several magnitudes. In the intervals between outbursts, 
which typically last years to decades, the X-ray luminos- 
ity is extremely low, and the optical flux is dominated by 
the luminosity of the secondary star (see Verbunt et al. 
1994, Van Paradijs & McClintock 1995 and Narayan, Mc- 
Clintock & Yi 1996 for reviews of the quiescence X-ray 
and optical properties of SXTs). Therefore, the properties 
of the secondary star can be studied during the quiescent 
episode. 

The SXT Nova Sco 1994 (GRO J1655-40), discovered 
on 1994 July 27 with BATSE (20-100 keV) on board th e 
Compton Gamma Ray Observatory (Zhang et al. 1994), 
has been studied extensively during the past two years 
in X-rays and at optical and radio wavelengths (Tin- 
gay et al. 1995, Bailyn et al. 1995a, Harmon et al. 1995b, 
Hjellming & Rupen 1995, Zhang et al. 1995, Bailyn et al. 
1995b, Paciesas et al. 1996, Van der Hooft et al. 1997). 
CCD photometry of Nova Sco 1994 revealed a double- 
waved m odulation of the o p tical light at a period o f 
2.6 days (|Bailyn et al. 1995bj |Van der Hooft et al. 19971) . 
Strong evidence that the compact object in Nova Sco 1994 
is a black hole was presented by Bailyn et al. (1995b) 
who established a spectroscopic period of 2.601 ± 0.027 
days and a mass function /(AT) — 3.16 ± 0.15 Mq. 
The secondary was cl assified as an F3-F6 IV type star 
( prosz fc Bailyn 1997p . 

As an SXT, Nova Sco 1994 is remarkable, due to 
the fast recurrence of its successive X-ray outbursts since 
the first X-ray detection of the source in 1994. The X- 
ray outbursts of Nova Sco 1994 were initially separated 
by intervals of about 120 days (Zhang et al. 1995) and 
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lasted many weeks (Harmon et al. 1995a, b; Wilson et al. 
1995). During these periods the optical brightness of Nova 
Sco 1994 was between 1 and 2 m agnitudes above it s 
quiescent V= 17.3 pre-outburst level ( Bailyn et al. 1995a ), 
which likely reflects the contribution of the accretion disk 
and X-ray heating of the secondary star. Most photomet- 
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ric observations of Nova Sco 1994 were made during this 
period. 

However, in the interval 1995 September to 1996 April, 
no significant X-ray flux was detected from Nova Sco 1994 
indicating a (temporary) longer period of quiescence. We 
present in this paper the results of extensive V, R, and 
i band photometry of Nova Sco 1994 obtained during 
1996 March, and derive its quiescence light- and color 
curves. We derive limits on the system parameters based 
on ~ 30 000 theoretical light- and color curves in the 
three photometric bands, calculated for a wide range of 
parameters of a model describing a Roche lobe filling sec- 
ondary star supplemented with a model of an accretion 
disk around the compact star. The limits on the mass 
of the secondary star are compared with theoretical evo- 
lutionary tracks which we computed for both single and 
binary star evolution. 



2. Observations and Data Reduction 

We observed Nova Sco 1994 (GRO J1655-40) during 28 
consecutive nights starting 1996 March 5 with the 0.91m 
Dutch telescope at the European Southern Observatory 
(ESO) in Chile. The telescope was equipped with a CCD 
camera and the observations were made using standard V 
and R filters and a Gunn i filter. Integration times were 
10 minutes for the V filter and 5 minutes for the R and i 
filters. We obtained 67 images of Nova Sco 1994 in the V 
band, and 69 in the R and i bands each. All images were 
corrected for the bias, and flat fielded in the standard way. 
A 3' x 3' i?-band find chart of the field of Nova Sco 1994 
(taken on 1996 March 5), in which Nova Sco 1994 and five 
nearby comparison stars are indicated, is shown in Fig. [l[ 

We applied aperture photometry to Nova Sco 1994 and 
five nearby comparison stars using MIDAS and addition- 
ally written software operating in the MIDAS environ- 
ment. The five comparison stars were checked for vari- 
ability over the entire data set. From the disperion of the 
relative magnitudes of each comparison star over the en- 
tire observing period (< 0.008 mag) we conclude that the 
brightness of these stars is constant. The accuracy of a 
single measurement of the differential magnitude of Nova 
Sco 1994 (relative to the average of the five comparison 
stars), taking into account Poisson noise only, is typically 
0.006 mag in each of the three bands. 

Photometric calibration of the data wa s performed us - 
ing images of two Landolt standard fields ( Landolt 1992 ), 
Rubin 149 (8 stars) and PC 1323 (4 stars), during nights of 
good photometric quality. Observations of different stan- 
dard fields within one night showed that the accuracy of 
the photometric calibration was better then 0.05 mag. 



3.1. Quiescence light and color curves of Nova Sco 1994 



During the period of our observations the average V mag- 
nitude of Nova Sco 1994 was 17.25 and the light curve 
showed a modulation with two equal maxima and un- 
equal minima, whose depths are 0.31 and 0.48 magni- 
tudes. The average V magnitude of Nova Sco 1994 is 
consistent with the pre-outburst brightness of V=17.3 
(Bailyn et al. 1995b), which indicates that the source was 
in quiescence during the period of our observations. Dur- 
ing two successive pointings on March 8, separated by ~ 5 
hours, we detected a very low optical brightness of Nova 
Sco 1994 in each of the three filters (for the F-band we 
then found V=17.48 and V=17.53, respectively). 

There is remarkably little variability superposed on the 
quiescence light curve of Nova Sco 1994 in comparison to 
several other quiescent SXTs, which show cycle-to-cycle 
variations of the order of a few percent (McClintock & 
Remillard 1986, Remillard et al. 1992, Martin et al. 1995, 
Chevalier & Ilovaisky 1996, Shahbaz et al. 1996). Haswell 
(1996) has recently suggested that such brightness varia- 
tions may be similar to those seen in supcrhumps during 
superoutbursts of SU UMa type cataclysmic variables. 

The data in the three bands were independently 
searched for periodicities using a Lomb-Scargle peri- 



3. Results 



odogram (Press & Rybicki 1989). The periodograms show 
five strong peaks, corresponding to half the orbital period 
and its one-day aliases, and associated with each peak the 
characteristic pattern resulting from the data sampling. In 
each band we found a distinct maximum at a fundamental 
period of 1.3123 days. Since the light curve of Nova Sco 
1994 features two minima with unequal depth, we folded 
the data on twice this period, using a trial zero point (cho- 
sen in the center of the data set). The time of inferior con- 
junction of the secondary, T; n f , was derived from the folded 
light curve, by fitting a parabola to the shallow minimum, 
which yielded a phase offset to be applied to the trial zero 
point. The shallow minimum occurs close to the epoch of 
inferior conjunction of the secondary star as predicted by 
the ephemerides of Nova Sco 1994 reported in previous 
investigations (see e.g. Bailyn et al. 1995b, Van der Hooft 
et al. 1997). 

The photometric period of 2.6246 days of Nova Sco 
1994 is slightly longer than previous measurements of the 
period based on both photometric and spectroscopic ob- 
servations (see e.g. Bailyn et al. 1995b, Van der Hooft 
et al. 1997). In order to improve the accuracy of the or- 
bital period we combined 9 photometric measurements of 
T; n f from our 1995 and 1996 data with photometric mea- 
surements by Bailyn et al. (1995b) and Orosz & Bailyn 
(1997), and two spectroscopic measurements by Casares 
(1996) and Orosz & Bailyn (1997). From a least-squares 
fit of this set of 13 arrival times to their respective cycle 
numbers (spanning 140 cycles) we derived the following re- 
fined ephemcris for the times of inferior conjunction of the 
secondary star: HJD 244 9838.4198(52) + 2.62168(14) x N, 
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with x^=2.3 for 11 degrees of freedom. In Fig. |^ we present 
the V, R and i light curves of Nova Sco 1994 folded at this 
ephemcris. The drawn lines represent light curves which 
were computed using a theoretical model of the secondary 
star which will be described in detail in Sect. 3.2. 

We calculated the (V-R), (R-i) and (V-i) color curves 
using consecutive V, R and i measurements taken close in 
time (< 20 min). These color indices, folded at the above 
ephemeris, are displayed in Fig. ||, together with theoreti- 
cal color curves (see Sect. 3.2). They show a clear modula- 
tion with the photometric period. The peak-to-peak am- 
plitudes in the (V-R) and (R-i) color curves are ~ 0.06 
and ~ 0.05 mag, respectively. The color curves have two 
unequal maxima and two equal minima. The maxima of 
the color curves occur near phase 0.0 and 0.5 (the largest 
occurs near phase 0.5), indicating that Nova Sco 1994 
is reddest near conjunctions; the system is bluest near 
quadratures (phases 0.25 and 0.75). 



3.2. Modelling the light curves: ellipsoidal model 

To construct theoretical light curves we calculated the 
flux of the secondary emitted into the direction of the 
observer as a function of orbital phase, by numerical inte- 
gration of the contributions from a large number of small 
elements on the secondary's surface. We used a numer- 
ical model, discussed in detail by Tjemkes et al. (1986), 
based on a description of the tidal and rotational distor- 
tion of the secondary in terms of an equipotential surface. 
The non-uniform surface brightness distribution across the 
secondary is described according to Von Zeipel's theo- 
rem. We adopted 0.25 as the value of the gravity dark- 
ening coefficient j3, valid for stars with a radiative enve- 
lope ( Von Zcipcl 1924 ). For each surface element we used 
a black-body approximation (at the local effective tem- 
perature) for the perpendicularly emerging specific inten- 
sity. We applied a standard linear limb darkening law, 
with limb darkening coefficients as a function of temper- 
ature taken from Al-Naimy (1978). For the theoretical 
spectra calculated by Kurucz (1979) one finds that over 
the effective temperature interval 5500 to 8000 K (cover- 
ing virtually all surface elements on the secondary) the 
average deviation of the gradient 9 log Fy/dT with re- 
spect to the blackbody approximation is ~ 6 % (similar 
results were obtained for the dependence of Fr and Fi on 
T [~ 3 — 4 %]). To estimate the effect of our blackbody 
approximation on our results we have made trial calcula- 
tions in which the Planck functions were multiplied by an 
arbitrary function of T which accounts for the difference 
with the Kurucz spectra. We find that the depth of the 
deepest minimum changes by Ss 2 % ( F-band) and 1 % 
[R and i-band) with respect to the blackbody approxima- 
tion, while the depth of the shallow minimum changes by 
£ 4 % ( ^-band) and <, 2 % (R and j-band). These differ- 
ences are small enough that their effect on the solutions 



in the (i, M%) plane (see below) can be ignored compared 
to the effects of uncertainty in the system parameters. 

Initially, the possible contribution of an accretion disk 
to the emergent flux was not included in our modelling, 
nor were mutual eclipses and the effect of X-ray heating on 
the secondary. In our calculations we used the secondary 
mass M2 and the orbital inclination i as free parame- 
ters. The observed mass function f(M) — 3.16 ±0.15 M© 
( Bailyn et al. 1995b| ) then fixes the primary mass Mi and 
the mass ratio q — M1/M2. Together with the bolometric 
luminosity L opt of the secondary, the temperature distri- 
bution across the latter's surface, and the light curve, are 
determined. 

We calculated a first set of 81 light curves (and their 
corresponding color curves) for i between 50° — 90° (steps 
of 5°), and secondary masses in the range 0.5 < M2 < 4.5 
M Q (steps of 0.5 M Q ). We adopted a bolometric lumi- 
nosity for Nova Sco 1994 of L opt =41 L Q , as determined 
from the mean ^brightness ( V = 17.25, see Fig. p|), a dis - 
tance to the system of 3.2 kpc (Hjcllming & Rupcn 1995), 
a bolometric correction B.C. = —0.03 mag (F5 V type 
star; Popper 1980) and a color excess E(B-V)=1.3 mag 
(Borne et al. 1996). Each of the theoretical light curves 



were given an offset such that the data observed at the 
shallow minimum (phase 0.0) was matched. We found that 
none of these light curves could reproduce the sharp and 
deep minimum at phase 0.5, indicating that our model re- 
quires additional structure. Near this phase the compact 
object is in front of the secondary star, and the light curve 
shows significant scatter whereas it is smooth during the 
remaining part of the orbital cycle. Both the sharpness 
of the minimum and the increased scatter near phase 0.5 
likely result from the presence of an accretion disk around 
the compact star, which is obscuring (part of) the light 
of the secondary at this phase. There is no evidence for 
a disturbance of the light curve near phase 0.0; we there- 
fore, in this first phase of our analysis, computed a xt 
for each set of parameters, using the light curves in the 
three photometric bands simultaneously from which the 
data obtained between photometric phase 0.35 and 0.65 
were excluded. 

We found that the only acceptable solutions are for i 
in the range 60° — 75°, for any of the secondary masses 
between 0.5-4.5 M©. We made a new set of theoretical 
light (and color) curves, using a denser grid of i and M2 
values: 61° < i < 74° (0?25 step size) and 0.2 < M 2 < 5.0 
M Q (0.1 M step size), respectively. 

From the fits of the theoretical to the observed light 
curves we find that the binary inclination and mass of the 
secondary star can be well constrained. The xt contours 
in the i versus M2 plane trace a narrow, ellipse-shaped re- 
gion. Such a correlation between i and M2 is expected as 
the amplitude of the brightness modulation will decrease 
with both smaller inclination angles and larger mass of 
the secondary star (the distortion of the secondary, and 
therefore the temperature contrast across the stellar sur- 
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face decreases as M2 increases). A minimal xl of 2.54 (130 
d.o.f.) was found at (67?25, 1.60 Mq), which indicates that 
minor variations are present in the light curve in the phase 
interval 0.65-1.35. The 1, 2 and 3<r confidence contours in 
the (i, M2) plane are shown in Fig. |4|, which were com- 
puted by scaling the required A% 2 by the minimum value 
of xl- 

The assumed value of the bolometric luminosity of the 
secondary star of Nova Sco 1994 (41 L Q ) is sensitive to er- 
rors in both the distance to the source and its color excess. 
The latter is probably the most important source of uncer- 
tainty as the distance is well determined at 3.2 ± 0.2 kpc 
( Hjcllming fc Rupcu 1995] ) based on a kinematic model of 



the radio jets; other distance estimates are consistent with 
this value [3.5 kpc (McKay & Kesteven 1994); 3-5 kpc 
(Tingay et al. 1995) both based on H 1 absorption; ~ 3 



kpc ( jBailyn ct al. 1995a ) based on interstellar absorption; 
- 3 kpc (|G reiner et al. 1995j ) based on a dust scattering 
halo] . A minimal value of the color excess was reported by 
Bailyn et al. (1995a) [E(B-V)=1.15 based on the equiv- 
alent width of the Na D-lines in optical spectra taken in 
1994 August], while a largest E(B-V) of 1.5 was reported 
by Home et al. (1996) [based on a power-law spectral fit 
to UBVRI photometry taken during 1996 May]; deep 220- 
nm absorption in a HST spectrum of Nova Sco 1994 take n 
in 1996 May suggested E(B-V) of 1.3 ( [Home et al. 1996|) . 
Therefore, we assumed E(B-V) to lie in the range 1.2- 
1.4, with the corresponding value of L opt in the range 31- 
54 Lq. We therefore also calculated theoretical light- and 
color curves for secondary luminosities of 31 (for which 
we took 59° < i < 74°), and 54 L (with 62° < i < 78°); 
the steps in i were 0?25. In all calculations we covered the 
range 0.2 < M2 < 5.0 Mq, respectively, with steps of 0.10 
M®. 

From these theoretical light curves we derived ellipse- 
shaped xl contours in the i versus M2 plane similar to 
those derived for the curves calculated for L opt — 41 L Q . 
However, the location of the contours in the (i, M2) plane 
depends on L opt . For 31 L Q the xl contours shift towards 
lower binary inclinations, whereas they shift to larger in- 
clination for 54 L . Minimal xl values (for 130 d.o.f.) were 
found at (63?75, 1.30 M ) [xl = 2.60; 31 L ], and (71?25, 
1.90 M Q ) [xl = 2.52; 54 L j. The 1, 2 and 3cr confidence 
contours were computed similar as in the case of 41 Lq, 
and are displayed in Fig. ||. 

3.3. Inclusion of an accretion disk to the model 

The model described in Sect. 3.2 quite well describes the 
data obtained at inferior conjunction of the secondary 
star, but not those obtained near superior conjunction, 
when we observe less light than the model predicts. A 
natural explanation for the flux deficit near phase 0.5 is 
that an accretion disk around the compact star obscures 
part of the secondary at certain viewing angles. In the 
second stage of our modelling efforts we included an ac- 



cretion disk in the model, described as a flat cylinder with 
a radial temperature distribution. The vertical outer edge 
of the disk is assumed not to emit radiation. By including 
such an accretion disk the number of parameters in the 
model is increased by three: the radius of the disk as a 
fraction of the effective Roche lobe, a, the flaring angle 7 
of the disk, and the temperature at the outer edge of the 
disk, Tcdgc- In view of the very low X-ray luminosity of 
Nova Sco 1994 during the period of our observations, we 
assumed that the radial temperature distribution in the 
disk follows an r~ 3 / 4 dependence (Pringle 1981). 

We computed ~ 25 000 theoretical light and color 
curves for intrinsic luminosities of 31, 41 and 54 Lq, a frac- 
tional disk radius between 0.6 and 1.0 (step size of 0.1), 
outer disk temperatures of 100 or 1000 K, and a flaring 
angle of the disk of either 2° or 10°. We used step sizes of 
0?25 and 0.1 M for i and M2, respectively. We performed 
several tests (using a less dense grid for i, and M2) in which 
temperatures at the outer edge of the disk were selected 
between 100 and 6000 K. Light curves computed with 
Todgc > 1000 K do not resemble the data well: for such 
outer edge temperatures the accretion disk contributes too 
much to the total flux of the system, which decreases the 
amplitude of the ellipsoidal modulation and changes the 
relative depth of the minima in the light curve. At supe- 
rior conjunction of the secondary star, the contribution of 
the disk to the total flux is < 4 % for T odgc = 1000 K, 
and is entirely negligible for lower values of T c dgc- Thus, 
the main effect of the disk is to occult a fraction of the 
secondary star near phase 0.5. 

In calculating xl we did not include the six deviat- 
ing points (two in each of the three passbands) obtained 
during subsequent pointings on March 8; they must be 
caused by variations in the disk structure, which cannot 
be incorporated into our model. 

For each set of parameters, the 1, 2 and 3cr confidence 
contours in the (i, Ma) diagram were computed as before. 
The confidence contours corresponding to models having 
an accretion disk with a flaring angle 7 = 2°, and an outer- 
edge temperature T cc i gc = 100 K, are shown in Fig. |[ The 
top-left panel of this figure displays the confidence con- 
tours for L op t = 41 Lq and fractional disk radii between 
0.6 and 1.0. These contours trace narrow (~ 1?5 wide) el- 
lipses in the (i,M2) plane. Since the light curve of Nova 
Sco 1994 is dominated by ellipsoidal variations (see Sect. 
3.2), the 'islands' derived for the different fractional disk 
radii are located in the same general area of the (i,M2) 
diagram as the contours which were derived from the pure 
ellipsoidal model (see Fig. |i|). As expected, for increasing 
values of a the confidence contours move toward smaller 
values of i. 

In the top-right panel of Fig. [|, the confidence con- 
tours for intrinsic luminosities of 31 (dashed), 41 (solid) 
and 54 L (dashed-dotted) are displayed, with the other 
parameters unchanged. From this panel it follows that for 
a given fractional disk size, and intrinsic luminosities be- 
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tween 31 and 54 the solutions move along a rather 
narrow track in the (i, M2) diagram. Solutions for larger 
secondary masses are found at smaller inclination angles, 
since an accretion disk of a given relative size can obscure 
the light of the secondary star at smaller values of i as 
M2, and therefore, the size of the secondary, increases. In 
the remainder of our analysis we have assumed that the 
fractional disk radius is in the range 0.7 to 0.9 (see e.g. 
Paczyhski 1977, and Frank et al. 1992). 

We can further constrain the system parameters by 
noting that the spectral type of the secondary limits the 
range of its effective temperature (averaged over the sec- 
ondary surface), and therefore of the secondary radius. 
Since the secondary is assumed to fill its Roche lobe (this 
fixes its average density) , for a given value of L opt this re- 
quirement limits the position of the system in the (i, M2) 
diagram. The effective temperature decreases with both 
smaller intrinsic luminosity and larger mass of the sec- 
ondary star (for a Roche lobe filling secondary the average 
density depends only on orbital period). The secondary of 
Nova Sco 1994 has be en classified as an F3-F6 type star 
flOrosz fc Bailyn 1997| ). According to Popper (1980), ef- 
fective temperatures of such stars are in the range 6330- 
6620 K. The M2 ranges corresponding to this range in av- 
erage effective temperatures have been indicated for each 
intrinsic luminosity individually in the top-right panel of 
Fig. m In the bottom-left panel of Fig. || this constraint is 
combined with those obtained from the fitting of the light 
curves. 

For L opt = 54 L Q , average effective temperatures of 
< 6600 K are obtained for secondary masses of 3.6 M Q and 
higher. This limiting mass is much larger than the maxi- 
mum mass of the secondary allowed by the disk models for 
this luminosity. Therefore, taking the mass limit imposed 
by the average effective temperature of the secondary star 
into account, we can exclude an intrinsic luminosity of 
54 L Q for Nova Sco 1994 (see top-right panel of Fig. |J). 
For the other two intrinsic luminosities, the average effec- 
tive temperature limits the secondary mass to the range 
1.6-2.1 M Q (31 L Q ) and 2.4-3.1 M Q (41 L Q ), well inside 
the collection of solutions for these disk models. The com- 
bined constraints for these two values of L pt have been 
indicated by hatched polygons in the lower-left panel of 
Fig. I 

Given the smooth variation of the location of the confi- 
dence contours in the (i, M2) plane in response to changes 
in the fractional disk radius or bolometric luminosity, we 
connected both hatched areas by smooth curves, thereby 
producing a final constraint on the system in the (i, Mg) 
plane, which is shown in the lower-right panel of Fig. |[ 
Note that the final constraint shown in Fig. || is derived 
for the case r c d gc = 100 K and 7 = 2° only. 

In Fig. H we have collected the results of similar anal- 
yses for four combinations of outer-edge temperature and 
flaring angle of the disk: (100 K; 2°) [top-left panel], 
(1000 K; 2°) [top-right panel], (100 K; 10°) [bottom-left 



panel], and (1000 K; 10°) [bottom-right panel]. This fig- 
ure shows that the solutions corresponding to models for 
Lodge = 1000 K (right panels) are not shifted much with 
respect to those for 100 K (left panels). The models with 
a flaring angle of 7 = 10° (bottom panels) are slightly 
shifted toward lower binary inclination with respect to the 
models for 7 = 2° (top panels). 

In Fig. and || theoretical light- and color curves are 
shown, computed for the parameters i — 68? 75, M2 = 
2.10 M , 41 L , a = 0.8, r o d g o = 100 K, 7 = 2°, typical 
for the collection of solutions we obtained. 



4. Discussion 

4-1. Quiescence light curves 

The average V magnitude of Nova Sco 1994 during our 
March 1996 observations was 17.25, consistent with the 
pre-outburst brightness of V=17.3 (Bailyn et al. 1995a), 



which indicates that the source was then in quiescence. 
This is confirmed by X-ray observations of Nova Sco 1994 
made with ASCA and ROSAT during March 1996, which 
showed that its X-ray luminosity then was 2 x 10 32 erg/s 
(Robinson 1997). No significant X-ray flux of Nova Sco 
1994 was detected with RXTE during the period of our 
observations (3a upper limit of 12 mCrab in the 2-12 keV 
band); X-rays were first detected from Nova Sco 1994 with 
RXTE on April 23, 1996 (Remillard et al. 1996). 

Also, the smoothness of the V, R and i light curves 
(spanning 11 consecutive cycles of Nova Sco 1994) is 
in contrast to earlier reports of significant cycle-to- 
cycle variations in the Nova Sco 1994 light curve dur- 



ing periods of known X-ray activity (Bailyn et al. 1995b, 
Van der Hooft et al. 1997). The relative depth of the two 
minima in the light curves of Nova Sco 1994 is consistent 
with a quiescence optical light curve caused by ellipsoidal 
variations, the deepest minimum occuring at superior con- 
junction of the secondary star as expected. The brighten- 
ing in X-rays in late April 1996 was followed by increased 
activity at optical, and radio wavelengths (Home et al. 
1996, Hunstead & Campbell- Wilson 1996, Hjellming & 
Rupen 1996). 

We determined a revised ephemeris for inferior con- 
junction of the secondary star with respect to the 
black hole by combining our measurements of T m { with 
those published in other investigations, spanning a to- 
tal of 140 orbits of Nova Sco 1994. This period dif- 
fers by 2.7cr from the period we derived previously 
(Van der Hooft et al. 1997), and 0.8ct from the spectro- 
scopic period derived by Orosz & Bailyn (1997). 

The colours of Nova Sco 1994 show a double- waved 
orbital modulation, with amplitudes of ~ 0.06 in (V-R), 
~ 0.05 in (R-i) and ~ 0.11 in (V- i). The system is red- 
dest near the conjunctions, bluest at the quadratures. This 
colour variation reflects the varying temperature across 
the surface of the secondary; its pronounced appearance 
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is the result of the wavelength and temperature depen- 
dence of the limb darkening, which gives relatively low 
weight to emission from near the apparent horizon of the 
secondary. 

The light curve of Nova Sco 1994 can be understood 
in terms of a model of a distorted secondary star in com- 
bination with an accretion disk. The disk eclipses part of 
the secondary near phase 0.5, but does not give a sig- 
nificant contribution to the total luminosity of the sys- 
tem. The largest variability in the light curve of Nova Sco 
1994 occurs close to phase 0.5, and is most likely caused 
by variability in the structure of the outer disk, causing 
variability in the fraction of the secondary that is being 
eclipsed. 

4-2. Mass determination 

Our modelling of the light curves has led to constraints 
on the system parameters, which have been summarized 
in Fig. ||, in the form of allowed regions in the (i, M 2 ) di- 
agram. These regions incorporate solutions over the full 
range in bolometric luminosity of the secondary. From 
Fig. ^ it appears that there is some dependence of the 
solutions on the assumed value of the flaring angle 7 of 
the disk, but that it is virtually independent of its tem- 
perature structure. The reason is that the light curve away 
from phase 0.5 does not allow the disk to contribute sub- 
stantially to the optical emission of the system, which sets 
an upper limit to its outer temperature of ~ 1000 K; for 
these temperatures the main effect of the disk is to eclipse 
the secondary. Since for a given value of i the fraction of 
the secondary that is eclipsed increases with 7, we find an 
anti-correlation between 7 and i for the allowed solutions 
(compare the top and bottom panels in Fig. ||) . 

Since the mass function of the system is relatively well 
determined, a region in the (i, M2) plane allowed by the 
light curves corresponds to a range of allowed values of the 
mass ratio and the mass of the compact object in Nova Sco 
1994. For 7 = 2° we find inclination angles in the range 
66?2 — 70?7 and secondary masses between 1.6 and 3.1 
Mq. This then implies a mass ratio, q, and mass of the 
compact object, Mi, in the ranges 2.43-3.99 and 6.29- 
7.60 M©, respectively. For 7 = 10° the allowed ranges of 
i and M 2 are 63?7 - 68? 1, and 1.6-2.45 Mq, respectively. 
The corresponding allowed ranges of q and Mi are 2.93- 
4.20 and 6.35-7.18 M , respectively. In Fig. | we present 
the region in the M2 versus Mi plane occupied by the 
collection of allowed solutions, for a flaring angle of the 
accretion disk of 2° (solid lines) and 10° (dashed lines). 

4-3. Theoretical evolutionary tracks 

Since the bolometric luminosity and effective temperature 
of the secondary are well constrained, it is possible to com- 
pare the position of the secondary star of Nova Sco 1994 
in the Herzsprung-Russell diagram (HRD) with theoreti- 



cal evolutionary tracks. We have calculated tracks for both 
single stars and stars in a binary system with a recent ver- 
sion of the evolution code developed by Eggleton (1971, 
1995). The equation of state and the opacity tables have 
recently been discussed by Pols et al. (1995). The mass 
loss in the binary calculations was assumed to be con- 
servative and in all calculations we used a metallicity of 
Z = 0.02. 

In Fig. [n] we present a HRD with theoretical evolu- 
tionary tracks of single stars with masses 2.0, 2.25, 2.5, 
3.0, and 4.0 M (solid lines). Based on this figure, we po- 
sition the secondary of Nova Sco 1994 on the track for a 
single star with a mass of ~ 2.2 Mq. Assuming a mass 
of 7 M Q for the black hole primary, we obtain an orbital 
period of about 2.6 days (by adopting a radius for the 
Roche lobe filling secondary from the evolution track at 
the center of the error box [see Fig. |l(J, applying Kepler's 
third law and the relation between orbital separation and 
effective Roche radius of a star). This value is in good 
agreement with the observed period of 2.62 days. How- 
ever, this is not the only evolutionary track that fits the 
observations. 

Since the secondary is losing mass, it must have been 
evolved from an initially more massive star. Therefore, we 
computed several evolution tracks of the secondary star 
in a black-hole binary system, for which we adopted in 
all cases a total system mass of 9.2 M Q . The effect of a 
different value for the total system mass (8-11 Mg) does 
not change the outcome of our results significantly. It fol- 
lows from Fig. [l^ that a secondary with an initial mass 
as large as 4.0 Mg can also evolve towards the current 
HRD position of the secondary of Nova Sco 1994. When 
the tracks of these stars intersect the error box bounded 
by 1.49 < logiopt/Lo < 1.73 and 3.80 < logT ofr < 3.82, 
the orbital period and mass are about equal to the case 
of an evolved single star. This is not surprising as the 
secondary is less massive than the accreting star, fills its 
Roche lobe and starts losing mass on the main-sequence. 
The mass loss then occurs on a nuclear time scale, as is 
the case for a secondary with an initial mass of 4 M Q . For 
the secondaries with lower initial masses, the mass loss 
starts further on the main-sequence or even at the moment 
when the secondary is already crossing the Hertzsprung- 
gap. Since the radius is expanding at a much higher rate, 
the time scale for mass loss is now smaller than the nu- 
clear time scale, but still a factor of 10 larger than the 
thermal time scale. For an initial mass of the secondary 

4.0 M Q , the present mass and orbital period for all 
tracks of the mass-losing star that cross the center of the 
error box shown in Fig. [l(] are between 2.0-2.3 Mq and 
2.6-3.1 days, respectively. Therefore, we conclude that the 
binary nature of Nova Sco 1994 is not very important for 
the determination of the present mass of the secondary 
based on evolutionary tracks. 

From standard evolutionary calculations for single 
stars, we find that the tracks of stars with a mass in the 
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range 2.1-2.4 M© are consistent with the current posi- 
tion of the secondary of Nova Sco 1994 in the HRD (see 
Fig. 0). This method to deter mine the mass of the sec- 
ondary star is entirely independent of the constraints de- 
rived in Sect. 3.3 from modelling the light curves of Nova 
Sco 1994. We emphasize that the mass range derived from 
evolutionary arguments is consistent with that obtained 
from modelling the light curves. 

4-4- Comparison with Orosz & Bailyn (1997) 

Our results are consistent with those obtained by Orosz & 
Bailyn (1997). These authors describe the quiescence light 
curve of Nova Sco 1994 (obtained during the same period 
as our data) in terms of an ellipsoidal model and an ac- 
cretion disk. The main differences between their analysis 
and the one discussed here are the value of the power-law 
exponent of the radial temperature distribution on the 
accretion disk (fixed to -0.75 in this work, free parameter 
in the analysis by Orosz & Bailyn 1997), and the use of 
the temperature at the pole of the secondary star as in- 
put parameter (fixed to 6500 K, Orosz & Bailyn 1997), in 
contrast to a range of intrinsic luminosities (31, 41, and 
54 Lq, this work) as input to our model. 

The best fitting (i, q) values of 69? 50 ± 0?08 and 
2.99 ± 0.08 (ler statistical errors) found by Orosz & Bai- 
lyn (1997), combined with the mass function imply a sec- 
ondary and black hole mass of 2.34 ± 0.12 and 7.02 ± 0.22 
Mq, consistent with our results. Also, the values for the 
fractional disk size (0.747 ± 0.010), and angular thickness 
of the disk (2?23±0?18) derived by Orosz & Bailyn (1997) 
are consistent with the results discussed here for the case 
of a thin (2°) accretion disk. However, our analysis shows 
that models with an accretion disk of 10° angular thick- 
ness are also able to describe the data well. The contribu- 
tion of the accretion disk to the total luminosity is gov- 
erned by its fractional size and mean temperature, and 
therefore, by the temperature at the outer edge of the 
disk and the radial temperature distribution on the disk. 
Orosz & Bailyn (1997) derived an outer edge tempera- 
ture of 4317 ± 75 K, significantly larger than the value 
we derived from our analysis. This probably reflects the 
different value of the power law exponent between both 
analyses; in our analysis the exponent was kept fixed to 
-0.75 (Pringle 1981), while Orosz & Bailyn (1997) treated 
the exponent as a free parameter and obtained a value of 
—0.12 ± 0.01. However, it should be noted that in both 
analyses the contribution of the disk to the total luminos- 
ity of the system is less than 5%, consistent with the value 
derived from spectroscopy during 1996 February (Orosz & 
Bailyn 1997). 

The errors quoted by Orosz & Bailyn (1997) are la 
internal statistical errors, derived from fitting the light 
curves of Nova Sco 1994. Note however, that in our anal- 
ysis we have included errors caused by systematic uncer- 
tainties in the properties of the system, e.g. in the inter- 



stellar absorption its light suffers. Therefore, we allow for 
a range of intrinsic luminosities of Nova Sco 1994 (31-54 
L Q ). We also allow a range of values for the parameters 
describing the size and temperature of the accretion disk 
(0.7 < a < 0.9,2° < 7 < 10°, and 100 < T eff < 1000 K), 
and adopt the 3a confidence contours (see Fig. ^j. There- 
fore, the limits to the system parameters we have derived 
in this study are somewhat larger compared to those ob- 
tained by Orosz & Bailyn (1997). 

4-5. Systematic effects 

In our modelling of the light curve of Nova Sco 1994 we 
have made a number of implicit assumptions, e.g., with 
respect to the geometry of the disk, and such assumptions 
may affect the results of our analysis systematically. There 
is one assumption whose systematic effect is particularly 
strong, i.e., that the non- uniform temperature distribu- 
tion across the surface of the secondary can be described 
by a standard gravity darkening coefficient (3 ~ 0.25. This 
value applies in a situation where the outer layers of a 
star are in radiative equilibrium. For stars with convective 
outer layers Lucy (1967) derived a gravity darkening coef- 
ficient of 0.08. At its spectral type (F3-F6) the secondary 
of Nova Sco 1994 is near the boundary in the Hertzsprung- 
Russell diagram separating hot stars with radiative en- 
velopes and cool stars with convective envelopes. There- 
fore, it may contain a shallow convective layer. We have 
investigated the applicability of standard gravity darken- 
ing by computing light curves for f3 in the range 0.08 to 
0.25 (in steps of 0.01). Decreasing (3 leads to a decrease 
in the brightness contrast over the secondary's surface, 
and therefore to a decrease of the amplitude of the light 
curve. Since the amplitude increases with i the strongest 
limits on (3 are obtained for i — 90°. For an assumed 
Mi = 2.2 M Q and L opt — 41 L Q (near the center of our 
solution space) we derive a lower limit (3 > 0.18 on the 
gravity darkening coefficient. This suggests that the as- 
sumption (3 = 0.25 is justified. 

5. Conclusion 

We conclude that during our March 1996 observations 
Nova Sco 1994 was in a state of true X-ray quiescence. 
We derive a refined ephemeris for inferior conjunction of 
the secondary star with respect to the compact star of 
HJD 244 9838.4198(52) + 2.62168(14) x N. We calculated 
theoretical light- and color curves, which constrain the bi- 
nary inclination and mass of the secondary star to lie in 
the range 63?7-70?7 and 1.60-3.10 Mq, respectively. The 
implied mass of the black-hole is in the range 6.29-7.60 
Mq. We find that an accretion disk is required to model 
the light- and color curve of Nova Sco 1994. The disk does 
not contribute significantly to the luminosity, and there- 
fore, acts as body of obscuration only. The mass range we 
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obtained for the secondary star is supported by the results 
of stellar evolution calculations. 
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Fig. 1. 3' x 3' _R-band find chart of the field of Nova Sco 1994, taken on 1996 March 5. The exposure time was 300 seconds; 
North is at the top, East to the right. Nova Sco 1994 is indicated by a circle and is numbered 1, the five comparison stars are 
numbered 2-6. 
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Fig. 2. V, R and i light curves of Nova Sco 1994 folded at the photometric 
ephemeris T inf (HJD) = 244 9838.4198(52) + 2.62168(14) x N. The curves are repeated over 1.5 cycles for clarity. Phase 0.0 
corresponds to inferior conjunction of the secondary with respect to the compact star. The drawn lines represent a theoretical 
light curve computed for a binary inclination and secondary mass of (67?25, 1.60 M ). 
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Fig. 3. (V-R), (R~i) an d (V-i) color curves folded at the photometric 

ephemeris T inf (HJD) = 244 9838.4198(52) + 2.62168(14) x N. The curves are repeated over 1.5 cycles for clarity. Phase 0.0 
corresponds to inferior conjunction of the secondary with respect to the compact star. The drawn lines represent a theoretical 
color curve computed for a binary inclination and secondary mass of (67? 25, 1.60 M ). 
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Fig. 4. 1,2,3a confidence contours in the i versus M2 plane for intrinsic luminosities of 31 (dashed), 41 (solid), and 54 L Q 
(dashed-dotted). The average effective temperatures from the secondary, calculated from the intrinsic luminosity and the effective 
radius of the secondary, are indicated for 6300 and 6600 K. 
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Fig. 5. 1, 2, 3a confidence contours in the i versus Mi plane for intrinsic luminosities of 41 [top-left panel], and 31 (dashed), 41 
(solid), and 54 Lq (dashed-dotted) [other panels], for accretion disk models with outer edge temperature T e d g e= 100 K, flaring 
angle 7 = 2°, and fractional disk radii between 0.6 and 1.0. Average effective temperatures of the secondary star are indicated 
in the top-right panel. In the bottom-left panel, the solutions are constrained by the smooth curves tangent to the a — 0.7, 0.9, 
and 31, 41 and 54 L Q confidence contours. The solutions in the i versus M2 plane which are constrained by these limits are 
indicated by the hatched area in the bottom- left panel. Finally, the total collection of solutions for this (T e d g e, 7) is given in 
the bottom-right panel. 
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Phase 

Fig. 7. Photometric data as presented in Fig. ^, from which six data points have been removed (see text). The drawn lines 
represent a theoretical light curve computed for a binary inclination and secondary mass of (68? 75, 2.10 Mq), an intrinsic 
luminosity of 41 Lq and a disk with a fractional size of 0.80, a flaring angle of 2° and a temperature at the outer edge of 100 K. 
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Fig. 8. Photometric data as presented in Fig. |3| from which six data points have been removed (see text). The drawn lines 
represent a theoretical color curve computed for a binary inclination and secondary mass of (68?75, 2.10 Mq), an intrinsic 
luminosity of 41 Lq and a disk with a fractional size of 0.80, a flaring angle of 2° and a temperature at the outer edge of 100 K. 
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Fig. 9. Collection of solutions in the M2 versus M\ plane for a flaring angle of the accretion disk of 2° (solid lines) and 10° 
(dashed lines). The mass of the secondary star is restricted to 1.6-3.1 M Q , while the mass of the compact object lies in the 
range 6.29-7.60 M©. 
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Fig. 10. Evolutionary tracks of single stars with mass 2.0, 2.25, 2.5, 3.0 and 4.0 M (solid lines) and for stars in binary system 
which lose mass to their companion (dashed/dotted curves). The total system mass of the binary system is in all cases 9.2 M . 



